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Abstract 



For the nonlocal T-periodic Gross-Pitaevsky operator, formal solutions of the Floquet problem 
asymptotic in small parameter h, h —* 0, up to 0(fi 3//2 ) have been constructed. The quasi-energy 
spectral series found correspond to the closed phase trajectories of the Hamilton-Ehrenfest system which 
are stable in the linear approximation. The monodromy operator of this equation has been constructed 
to within 0(h 3 ^ 2 ) in the class of trajectory-concentrated functions. The Aharonov-Anandan phases have 
, been calculated for the quasi-energy states. 

o 

CN ■ 1 Introduction 



The global properties of the solutions of mathematical physics equations which simulate a physical system 
are determined by the nontrivial geometry and topology of the system. A basic problem in investigating the 
■ behavior of a system is to find mathematical structures which would describe the topology of the system 
and make possible its efficient analysis. For quantum systems subjected to an external action cyclically 
varying with time, this structures are well known as topological or geometric phases (GP's) of the wave 
function. Some manifestations of GP's were also observed in polarization optics and mechanics. Theoretical 
and experimental studies of geometric phases in quantum mechanics have been performed since the late 20s 
of the last century. The work by Berry [1] has substantially extended the area of application of GP's and 
has led to the interpretation of this notion in terms of gauge symmetry and an effective gauge field or, in 
geometric formulation, in terms of a Hilbert bundle with a finite-dimensional base. A detailed description 
jJJ ■ of this problem can be found, for instance, in [2-4]. 

The theory of GP's in quantum mechanics is based on the linearity of the Schrodinger equation. For 
nonlinear equations, the notion of GP can also be introduced (see, e.g., [2]), but the GP's in nonlinear 
systems are less well understood. Difficulties arise not only for lack of the principle of superposition of 
solutions. In case of a nonlinear equation, the expression for the GF's is constructed by analogy with the 
linear theory (see, e.g., [2]), and it is not fully obvious that the expression constructed is governed only by 
the geometry of the system and that it does not involve the dynamic contribution due to nonlinearity. The 
nontrivial topology of a system can be determined by the boundary conditions as well as by the external 
fields. The latter appear in the equation as variavle factors. Few, if any methods are available to construct 
exact solutions for equations of this type. Therefore, it is natural to study GP's in nonlinear systems by 
using adequate approximation methods. 

In this connection, the analysis of the GP's for the nonlocal Gross-Pitaevsky equation given below 
would be of interest. The evolution of the initial state of the nonlinear equation in the semiclassical 
approximation is determined in fact by a set of associated linear Schrodinger equations. This allows one not 
only to calculate in full measure the GP for the nonlinear equation in explicit form, but also to interprete it 
in terms of "an induced gauge field" (see [2]), which in turn proves its geometric origin. 

A class of solutions asymptotic in small parameter h — > 0, which are localized in a neighborhood of some 
phase curve, has been constructed [5-9] for the nonlocal Gross-Pitaevsky equation (NGPE) with an external 
field (variable coefficients) and a nonlocal nonlineariry 

{-ih8 t +'fl(t)+HV(jL,V)}9 = 0, $6L 2 (K"). (1.1) 
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Here, the operators 

H(t) = H(z,t), (1.2) 

V(t,*) = J dy**(y,t)V(z,w,t)V(y,t) (1.3) 

are functions of the noncommuting and Weyl-ordered operators 

the function 'J* is complex conjugate to f , x is a real parameter, and h is a "small parameter", h E [0, 1[. 
For the operators z and w the following commutation relations are valid: 



[zk,Zj]- = [wk,u>j]- = ihJ k j, [zk,Wj]- = 0, k,j = l,2n, (1.4) 
where J = \\JkjW2nx2n is a simplectic unit matrix 
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The goal of this work is to construct explicit expressions for the quasi-energy spectral series corresponding 
to the closed phase trajectories of the Hamilton-Ehrenfest system which are stable in the linear approxima- 
tion and GP's in the class of asymptotic solutions of equation (jl.ip - (jl.3p . constructed in [5-9], for the case 
where the operators Tt(t) and V(t, 1 J r ) are T-periodic functions of time: 

H(t + T)=H{t), V(t + T,V) = V(t,V). (1.5) 

For the linear Schrodinger equation(>f = 0) with a T-periodic Hamiltonian H(t), Zeldovich [10] and 
Ritus [11] were first to introduce an important class of solutions — quasi-energy states ^e(x, t, h), which 
possess the property 

^ £ {x,t + T,h) = e- i£T ^ £ {x,t,h). (1.6) 

The quantity £, involved in (|1.6p . is called quasi-energy, and it is defined modulo tko (w = 2ir/T), i.e. 
£' = £ + mfhu, m £ Z. States of this type play the key part in describing quantum-mechanical systems 
under strong periodic external actions for which the conventional methods of nonstationary perturbation 
theory are not applicable. 

For the nonlocal Gross-Pitaevsky equation (jl.l[) - (|1.3[) . we shall set the problem of finding a quasi-energy 
spectrum in the form (|1.6j) . 

The quasi-energy states (|1.6|) are a particular case of the cyclic states introduced by Aharonov and 
Anandan [12] (see also [2]). By cyclic evolution of a quantum system on a time interval [0, T] is implied the 
following: the state vector \I/(i) has the form 

(i) = e i/(t V(*), ie[0,T], (1.7) 

where 

f(T) — /(0) = <fi(mod 2n), <p(T)=<p(0). (1.8) 
The total phase (f> of the function al. 7[) is subdivided into two terms: the dynamic phase 

d --nJ dt mm) (L9) 



and the Aharonov- Anandan geometric phase 

7 = ,LMMf. (i.io) 


Comparing (|1.6p and Q1.7p . we obtain that for the quasi-energy states the function f(t) is given by 

/(*) = -£t/h, (1.11) 



while for the total phase (f>, according to (|1.8p , we have 



£T 



(mod 2ir). 



(1.12) 



By virtue of (|1.9p ~ (|1.12p , the Aharonov-Anandan phase 7^ corresponding to a given quasi-energy state 
^£(x,t,H) can be determined by formula (mod 2tt): 



-JJL + \ f ; A*e\Mt) + xV(tMt))]\*e) 
l£ h + hj (Hf e \9 e ) 




(1.13) 



2 Statement of the problem in the class of semiclassically 
concentrated functions 

We shall construct asymptotic solutions for the nonlocal Gross-Pitaevsky equation with the following 
assumption for the Weyl symbols of the operators Tt(z,t) and V(z,w,t) in (|1.1[) - (|1.3|) : 

Assumption 1 The functions 7i(z,t) and V(z,W,t) are C°° -smooth functions and for any multiindices a 
and (3, (a,j3 € Zfj_ n ) and T > 0, there exist constants k > 0, C a (T), and C a p(T), such that 



d^H(z,t) 



dz c 



<C a (T)(l + \z\r, 



d^ a+ ^V(z,w,t) 



dz a dw< 3 



<C Q/3 (T)(l + |z| + M) K , z,w€ 



Here. 



< t < T. 



d^V(z) 
dz* 



dz^dz* 2 ...dz2 n 



u,j G Z +! j = 1, 2n, 



(ai, a 2 , ■ ■ ■ , a 2n ), |a|=ai+a 2 



a2r, 



We now turn to the description of the class of functions in which we shall seek localized asymptotic 
solutions of the equation (jl.l|) - (|1.3|) . The functions of this class, which singularly depend on a small parameter 
h, represent a generalized notion of a solitary wave. They depend on an arbitrary phase trajectory z = 

z(t,h) eK£», t g 

with k = 0). In the limit h 
moving along a given phase curve z = Z(t, 0). Functions of this type are well known in quantum mechanics. 
In particular, among these are the coherent and "squeezed" states of quantum systems with a quadratic 
Hamiltonian [13,14]. 

Let us denote this class of functions by V^(Z(t, H),S(t, h)) and define it as follows: 



1 and on a real-valued function S(t, h) (analog of a classical action for a linear case 
0, the functions of this class are concentrated in a neighborhood of a point 



Vl = Vi{Z(t,K),S(t,H)) = 

= |$ : $(£ , t, h) = Vhj exp ^(S(t, K) + (P(t, ft), Ax)) 



(2.1) 



where the complex- valued function (p(£,t,VK) belongs to the Schwarz space S in the variable £ G K™, 
smoothly depends on t and regularly depends on \f~h as li — > 0. Here, Ax = x — X(t,h), and the real 
function S(t,h) and 2n-dimensional vector function Z(t,H) = (P(t,H),X(t,H)), which characterize the class 
Vfi{Z(t,h),S(t,h)), regularly depend on \fh in a neighborhood of h = 0, such that S(0,h) — 0, Z(0,h) — 
zq = (po) xq) is an arbitrary point of the phase space Rp™, and are to be determined for t > 0. 

The functions belonging to the class P^, at any fixed point in time t G M 1 , are concentrated, in the limit 
h — > 0, in a neighborhood of a point lying on the phase curve z = Z(t,0), ( e I 1 (the exact sense of this 
property is specified below by formulas ()2.6|) . (|2.7|) ). Therefore, it is natural to give the functions of the class 
P l h the name trajectory- concentrated functions as h — > 0. The definition of the class of trajectory-concentrated 
functions involves a phase trajectory Z(t, H) and a scalar function S(t, K) as "free" parameters. It turned 
out that these "parameters" can be uniquely determined by the Hamilton- Ehrenfcst system (see Sect. 3), 
which corresponds to the nonlinear (x 7^ 0) Hamiltonian of equation (|1.1[) - (|1.3|) . Note that for a Schrodinger 
type equation, in the linear case (h = 0), the vector function Z(t,0) — the principal term of the expansion 
in h — > — determines the phase trajectory of a Hamilton system with a classical Hamiltonian TL{p,x,t), 
and the function S(t,0) is a classical action along this trajectory. In particular, in this case 
involves the well-known dynamic (squeezed) states of quantum systems with quadratic Hamiltonians 

Let us consider the basic properties of functions of the class P l n {Z(t, K), S(t, h)) (see [5,6, 15]). 



the class P l h 



1 For functions of the class "Pj^Z (t, H), S(t,h)), the following asymptotic estimates are valid for the centered 
moments A a (t, H) of order \a\, a £ Z^™: 



A a (t,H) = A*(t,h)= =0(hW 2 ), h^O, H^O; (2.2) 

A a (t,h) = i, |a| = o. 

Here {Az} a denotes an operator with the Weyl symbol (Az) a , 

Az = z — Z(t, h) = (Ap, Ax), Ap=p- P(t,h), Ax = x-X(t,h). 

Denote by the symbol <J(h k ), k > 0, an operator F, such that for any function <E> £ V\{Z(t,K), S(t,H)) 
the following asymptotic estimate is valid: 

\\F$\\ i 

= o(h k ), h^O. (2.3) 



2 For functions of the class Vt(Z(t,K), S(t,K)), the following asymptotic estimates are valid: 

{-ihd t - S(t,h) + (P(t,h),x(t,h)) + (Z(t,h),JAz)} = 0(h), (2.4) 
{Az} a = 0(ftH/ 2 ), aeZf, h^O, 

and, in particular, 

Ax k = D(Vh), Apj = 0(Vh), k,j = T~n. (2.5) 

The following property gives an exact sense to the notion of the concentration, as h — ► 0, in a neighborhood 
of a point on the phase trajectory for functions of the class V\- 

3 For any function <fr(x,t, K) £ V\(Z(t, h), S(t, H)), the following limiting relations are valid: 

lim pp|$(f, t, K)\ 2 = S(x - X(t, 0)), (2.6) 
Iim^l^, t, h)\ 2 = S(f- P(t, 0)), (2.7) 
where Q(p,t,h) = F%^^g^(x,t,K), F^g_ >p is a Fourier h-transform [16]. 

4 Denote by (L(t)) the mean value of an operator L(t), tel 1 selfadjoint in L2(K") that is calculated by 
the function &(x, t, h) £V\ l . Then, for any function &(x, t, h) £ V\(Z(t, K), S(t, H)) and any operator A(t, h) 
whose Weyl symbol A(z,t,H) satisfies Assumption^ the following equality is valid: 

lim(A(t,h)) = lim -^($(x,t,h)\A(t,H)\<f>(x,t,H)) = A(Z(t, 0), t, 0). (2.8) 
ft— >o n^a <P 

By analogy with the linear theory [15], we give the following definition. 

Definition 1 A solution $>(x, t, K) of equation (II. \\ - (|1.3|) in the class of functions V l n is called semiclassically 
concentrated as h — > on the phase trajectory Z(t, h). 



The limiting character of the conditions (|2.6[) . (|2.7[) and the asymptotic character of the estimates 
(12. 5[) . valid in the class of trajectory-concentrated functions, make it possible to construct the semiclassically 
concentrated solutions of the nonlocal Gross-Pitaevsky equation not exactly, but approximately. In this 
case, the L2-norm of the error is of order h a , a > 1, as h — > on any finite time interval [0, T]. 

Denote such an approximate solution by 'fas = ^^(x, t, h). It satisfies the following problem: 

+ #(*) + *V(t, *as)l *as = 0(H a ), (2.9) 

at J 
* as £V t n (Z(t,h),S(t,h),h), t£[0,T], (2.10) 

where 0(h a ) denotes the function g^ a '(x,t,h) with represents the error for equation (|1.1[) ~ (|1.3[) and obeys 
the estimate in the norm L^: 

ms^ T \\g i - a \x,t,h)\\=0(h a ), H^O. (2.11) 



We shall also call the function v E' as (a;, t, K), satisfying the problem (|2.14p - (|2.1ip , a semiclassically concentrated 
solution (modh a , h — > 0) of the nonlocal Gross-Pitaevsky equation (|1.1[1 - (|1.3[1 . 

Thus, the semiclassically concentrated solutions 9 m (oS, t, K) of the nonlocal Gross-Pitaevsky equation 
approximately describe the evolution of the initial state ^(x, h) if it is chosen in the class of the trajectory- 
concentrated functions V\: 

9{g,t,h)\t=a = M^ ^06^(20,0). (2.12) 
The functions of the class V\ have the form 

^,h)=exp^{^o,x-x )}}ipo(^^,Vh), ^o&v'fiJeS^), (2.13) 

where zq = (po, xq) is an arbitrary point of the phase space M 2 ,". 

As in the linear case [17], among the solutions of equation (jl.ip that satisfy the quasi-periodicity condition 
(II. 6p . in the class of solutions concentrated in a neighborhood of some phase trajectory, we can isolate a set 
of semiclassical asymptotics which possess the following properties: 

1) ^ £ u (x,t,K) are approximate in mod 0(h 5 / 2 ) solutions of equation JO}. This implies that 

[ - ihd t + Hit) + xV(t,^(t))] (x, t, K) = v v (x, t, h), 

max \\v v (x,tM\L^ q )=0(h*l 2 ). (2J4) 

2) The functions "i>£ u have the form of wave packets concentrated at any time t in a neighborhood of a 
given T-periodic trajectory z = Z(t, H). We shall call these states, if any, the quasi-energy trajectory-coherent 
states (TCS's) of the nonlocal Gross-Pitaevsky equation. 



3 System of Hamilton— Ehrenf est equations 

The symbols H(z, t) and V(z, w, t) satisfy the conditions of Assumption[TJ Therefore, the operator H(z,t) 
in ljl.l|) - (|1.3j) is selfadjoint with respect to the scalar product in the space L 2 (M.™) and the operator 

V(z,w,t) with respect to the scalar product L 2 (M.^): (*(t)|$(t)) K 2n = / dxdy$*(x,y,t,h)$>(x,y,t,K). 

Hence, the squared norm of the exact solutions of equation (jl.ip - (ll.3p is reserved: ^(i)!) 2 — ^(O)!) 2 , and 
for the mean values of the operator A(t) = A(z, i) calculated for these solutions 

±(A(t))=(^)+ l -([H,A(m (3.1) 

where [A, B] = AB — BA is the commutator of the operators A and B. By analogy with a linear (x = 0) 
Schrodinger equation in quantum mechanics, we call equality (|3.ip the Ehrenf est equation for the mean values 
of the operator A(t) which corresponds to the nonlocal Gross-Pitaevsky equation. 

We suppose that for the nonlocal Gross-Pitaevsky equation (|l.ip - (|1.3p there exist exact (or differing by 
0(h°°) from exact) solutions in the class of trajectory-concentrated functions. Let us write the Ehrenfest 
equations (|3.ip for the mean values of the operators z, {Az} a , which are calculated by this type (trajectory- 
concentrated) solutions of equation (|l.ip - (|1.3p . with the use of the composition rules for the Weyl symbols 
[16]: 

where C(z) is the symbol of the operator C — AB and the numbers 1 and 2 above an operator indicate 
the order of its action (recall that z = (p,x), Z(t,H) = (P(t,K),X(t,K)), Az — z — Z(t,H)). Then, after 
calculations similar to that in the linear case x = (see for details [15,18]), restricting ourselves to the 
second-order moments, we obtain the following system of ordinary differential equations: 

i = Jdz (l + \(dz, A 2 a 2 ) + ^(d w ,A 2 d w )) (H(z, t) + xV(z, w, t)\ w=g ), 
A 2 = Jfi zz (z, t)A 2 - A 2 Sj zz (z, t)J, 
where A 2 — (A^ 1012 ) is a (2n x 2ra) matrix of "variances", the (2n x 2n) matrix fi zz (z, t) is defined by 



!S„(z,t) = [H zz (z,t) + xV zz (z,w,t)\ ]. 



and x — x\\ipo(x, ft)|| 2 , ipo(x, K) denotes the initial condition of the Cauchy problem (|2.12[) . 
System (I3.2p can be written in an equivalent, setting in the second equation 

A 2 (t)=A(t)A 2 (Q)A + (t), 

which then takes the form 

A = JSj zz (z,t)A A(0)=I. (3.3) 

We shall call the system of equations Q3.2p . (|3 . 3|) the second order Hamilton-Ehrenfest (HE) system (here 
M = 2 is the order of the greatest moment taken into account) corresponding to the Gross-Pitaevsky 
equation(fTl ^ -(fL3] ) . 

Theorem 1 Let functions ^^ 2 \x,t,h) and $ (x , t , K) be, respectively, an asymptotic, to within 0(h 3 ^ 2 ), and 
an exact solution of equation (|2.1|) . which coincide at some point in time t = to- Then the second order 
Hamilton-Ehrenfest systems (|3.2p . constructed for the functions ^( 2 \x,t,K) and ^(x,t,h), coincide. 

Theorem 2 Let functions ^^ 2 \x,t,K) and^{x,t,K) be, respectively, an asymptotic, to within Olh 3 / 2 ), and 
an exact solution of equation (|2.ip , which coincide at some point in time t = to . Then 

A 2 *(t, h) = A 2 ^ (2) (t, h) + 0{h 3 ' 2 ) = A 2 (t, h) + o(h 3 / 2 ), 
z 9 (t,h) = z^ (2) (t,h) + o(h^ 2 ) = z 2 (t,h) + o{h 3 ' 2 ). 



Let us consider the HE system (|3.2p as an abstract system of ordinary differential equations with arbi- 
trary initial conditions. Obviously, not all solutions of the HE system (|3.2p can be obtained by averaging 
the corresponding operators over the solutions of the Gross-Pitaevsky equation (jl.ip . For instance, the 
mean values should satisfy the Schrodinger-Robertson uncertainty relation [19], which, for the second order 
moments implies nonnegative determinacy of the matrix 

ih 

A 2 (t)--J. (3.5) 

In the one-dimensional case, condition (|3.5p is equivalent to the well-known Schrodinger uncertainty relation 
[20] 

h 2 

VppVxx - crl p > —. (3.6) 
Denote by $(t, £) the general solution of the system of Hamilton-Ehrenfest equations (|3.2p : 

g(t, C) = (P(t, h, q,X(t, h, £), A u (t, h, £), A 12 (t, «,£),..., A 2n2n (t, h, <£)) T (3.7) 
and by q the operator column 

0= (p>,x,{Ap 1 ) 2 ,Ap 1 Ap 2 ,...,{Ax n ) 2 ) 1 . (3.8) 

Here, 

£ = (Ci,...,CV) T el 3n+2 " ! (3.9) 

are arbitrary constants, and B 1 denotes the transpose to the matrix B. 

These constants £, Eq. (|3.7p specify the trajectory of a point in the phase space Ai N . 

Lemma 1 Let ^(x, t) be a particular solution of equation (|2.ip with the initial condition ^(x, t)\ t _ Q = ip(x). 
Determine the constants <t(^(t)) from the system 

fl (t, C) = J^(*(t)m®) + 0(h 3 / 2 ), (3.10) 

and the constants £("0) from the system 

8(o, c) = o(n s/2 ). (3.ii) 

Then 

£(¥(«)) = £(<ip) + 0(h 3/2 ); (3.12) 
that is, <t(^(t)) are asymptotic, to within (0{h 3 / 2 )), integrals of motion for equation (|2.ip . 



Proof. By construction, the vector 



1*1 



(*(t)|fl|*(i)) = fl (t, £(*(*))) 



(3.13) 



is a particular solution (to within 0{h z / 2 )) of system (|3.2p . which coincides with g(£, (£(■!/>)) at £ = 0. By 
virtue of the uniqueness of the Cauchy problem for system (|3.2[) . the relation 



(£, £(V0)=fl(t, £(*(*))) + 0(ft 3 / 2 ) 



(3.14) 



holds true. 

By virtue of the estimates (|2.2p , the Hamilton-Ehrenfest system can be solved not exactly, but approxi- 
mately: 

Z(t) = Z (t) + hZ x {t) + o(h 2 ), 

(3.15) 

A 2 (£) = A 2 Q) (t)+0(h 2 ). 

Substituting (|3.15[) to system (|3.2p . we obtain, to within 0(£i 3 / 2 ), a system of equation for zq — Zq, z\ = Z\, 



A 2 = A< 0) (£) 



z = Jd ZQ [H(z Q ,t) + >cV(z ,w,t)} 
z\ = J[Sj zz (z ,t) + kV zw (z a ,z ,t)]z 1 + 

+ — Jd z „Sp{[H zz (z ,t) + xV zz {z Q ,w,t) + xV ww (z ,w,t)]A 2 } 
[A 2 = JSj zz {z ,t)A 2 - A 2 Sj zz (z ,t)J. 



(3.16) 



The first equation of (|3.16p is a generalization of the system of Hamilton equation to the case of a self-action. 

If a solution zq = Zq (£) of this system is known and also known is the set of solutions of two systems in 
variations (for this to be the case, it is sufficient that the matrix V zw (z, z, £) be symmetric) 



a k = JS) zz (t)a k , k = 1,71, S) zz (t) = Sj zz (Z (t) 7 t) 
normalized by the condition 

{a fe (£),a,(£)} = K(£),<(£)} = 0, {a k (t), <(£)} = 2iS kl , 
where {a, b} is a skew-scalar product: 

{a, b} = {J a, b) = (W a , Z b ) - (Z a , W b ), 



(3.17) 
(3.18) 
(3.19) 



W a 
Z a 



Z b 



and of the system 

h k = Jft zz (t)a k , k = T^, fi zz (t)=f) zz (Z (t),t) + kV zw {Z (t),Z (t),t), 
also normalized by (|3.18|) . Then the general solution of the two last equations of (|3 . 1 6[) has the form 

n 

= £[&*(*)<»*(*) + 63E(*K(*)]. 

k=l 



A 2 (t) = A(t)VA^(t), 



(3.20) 

(3.21) 
(3.22) 



where 



o 

F(t) = ^Jd z Sp{ [H zz (z,t) + xV zz (z, w, £) + kV ww (z, w, £)] A 2 } 
A(t) = ( ai (£), a 2 (£), . . . , a„(£), aj(t), a* 2 {t), <(£)). 



w—z—Z{)(t) 



(3.23) 



Here, B k are integration constants and V is an arbitrary constant matrix. Thus, in this approximation, the 
solution of the problem is completely determined by the solution of the generalized Hamilton system and 
the system in variations. 



The quasi-periodicity condition (|I.6p leads to the following constraint on the solutions of the Hamilton - 
Ehrenfest system: 

Z (t + T) = Z (t), Z 1 (t + T) = Z 1 (t), A 2 (i + T)=A 2 (t) (3.24) 

We shall denote by Ct the values of the constants £ at which condition (13.24p is fulfilled. 

In this case, the systems in variations (|3.17p . (|3.20p are systems of ordinary differential equations with 
periodic coefficients 

Si xx (t + T) = Si„(t), f) zz {t + T) = fl„(t). 
State for the systems in variations (|3.I7p . (|3.20p the Floquet problems: 

a k (t + T) =e in « T a k (t), a k (t + T) = e^« T a k (t). (3.25) 

Suppose that each Floquet problem of (|3.17p . (I3.20p . (I3.25P has n, n = dimM™, linearly independent 

Floquet solutions a k (t) = (W k (t), Z k (t)) J , a k (t) — (W k (t), Z k (t)y , with pure imaginary Floquet indices 
iQ n ,iQ n , which satisfy the orthogonality and normalization condition (|3.18[) . 

It should be stressed that in the context of the Floquet theory for linear Hamilton systems with periodic 
coefficients, the conditions Im Q k ~ imply stability of the phase trajectory z — Z (t) in the linear 
approximation. Recall that the 2n-dimensional vectors a k (t) and a£(i), k = l,n, constitute a simplectic 
basis in C„", and the n-dimensional plane r n (Z(t, h)) spanned over the vectors a k (t) constitutes a complex 
germ on Z(t,h) [21,22]. 

Determine the integration constants B k from the condition of periodicity of the function hZ\(t) in time: 

Zx{t + T) = Z x {t). (3.26) 

Like for (|3.2ip . we put 

t 

h(t) = | J dr{F(r), o2(r)} + B k . (3.27) 
o 

Then conditions (j3~2"5)) , (|3~2fJ]l imply that 

b k (t + T) = b k (t)e~^ T . (3.28) 

Note that 

t+T t+T 

h(t + T) = Yi J dT{F(T\al(T)} + B k = b k (T) + ±: J dT{F( T ), a * k (T)} = 

T 
t 

= b k (T) + 1 J dr{F(r + T), 4(r + T)} = b k (T) + er^ T (b k (t) - B k ). 
o 

Hence, 

B k =b k {Ty^ T . (3.29) 



4 The linear associated Schrodinger equation 

The linearization of the nonlocal Gross-Pitaevsky equation in the class of trajectory-concentrated func- 
tions is the central point in our approach. 
Introduce the notation (see l|2.2p ) 

fl*(*) = ]f^jp<*(*)lfll*(*)) 

and expand the "kernel" of the operator V(t, \I r ) in a Taylor series in powers of the operators Aw = w — 
z^(t,h). Substituting this expansion in equations (| 1 . 1 p - (|I.3D . for the functions € V\, in view of the 
estimate (|2.2I) . we obtain 

LW(t,V)y = \-ihd t + Sj{z,t) + 3cV(z,w,t) + ^Sp[V ww (z,w,t)A 2 *(t,h)]}\ * = 6(h 3 / 2 ). (4.2) 



Let us associate the nonlinear equation (|4.2p with the linear equation that is obtained from (|4.2p by 
substituting the corresponding solutions of the HE system g(t, €t) for 0*(t), the mean values of the operators 
of coordinates, momenta, and second order centered moments. As a result, we obtain the following equation: 



L (2) (*,£t)$ = 0, $€P£, 



L {2) (t, € T ) = -ihd t + H(z, t) + xV(z, w, t) 



W = Z{t,<tT) 



V ww (z, w,t) 



w=Z(t,€ T ) 



A 2 (*,Ct) 



(4.3) 
(4.4) 



By virtue of (|3.24p . equations (|4.4p are equations with periodic coefficients. The quasi-periodicity condition 
(|1.6[l reduces, for equation (|4.4p . the Floquet problem 



$s(x,t + T, ft) 



_ e -iST/h 



<f>e(x,t,h). 



(4.5) 



Thus, the change of the quantum means of the operators w and {Aw} a by the solutions of the second 
order Hamilton-Ehrenfest system in equation (14. 2p linearizes the nonlocal Gross-Pitaevsky equation Ijl.ljl - 
()1.3p to within 0(ft 3 / 2 ). Hence, to construct the semiclassically concentrated states mod ft 3 / 2 of the nonlocal 
Gross-Pitaevsky equation Ql.l|) - (|1.3|) . it suffices to construct, with the same accuracy, an asymptotics of the 
solution of a Schrodinger type linear equation. 



Definition 2 We call equations of the form (14. 3p parametrized by constants £ (|2.12p the set of associated 
linear Schrodinger equations for the nonlocal Gross-Pitaevsky equation (jl.l|) - (|1.3|) . and denote by <£> = 
$(x, t, H, £) its (arbitrary) solution in the class of functions Vt Q 



The following statement is obvious. 



Statement 1 If a function $>^(x,t, h, £) 
equation (|4. 3[) — satisfies the initial condition (|4.5 



symptotic (to within 0(ft 3 / 2 ), h — > 0) solution of 



^ 2 \x 1 t 1 h,^ ))\ t=0 = iJo, 



(4.6) 



£ften £fte function ^^(x,t : h) = <I>( 2 ) (x, t, ft, £(-0o)) is o asymptotic (to within 0(ft 3 / 2 ), ft — > 0) solution of 
the Cauchy problem for the nonlocal Gross-Pitaevsky equation (jl.ll) . Constants €(ipo) are determined by the 
equation (|3.1ip . 

To construct the function <I>( 2 )(x, £, ft, -0o) G ^(-^(^ we shall use the technique for construct- 

ing localized asymptotics developed in [17, 21-24] for linear quantum-mechanics equations, modifying this 
technique taking into account the initial statement of the problem (|2.14p - (|2.12p and the structure of the 
Hamiltonian (|4.4p for the (linearized) nonlocal Gross-Pitaevsky equation in the trajectory-coherent approx- 
imation. Namely, we shall define the point of localization in M. 2n x of the semiclassically concentrated solution 
of equation (|4.3p — "parameter" Z(t, ft) entering into the definition of the class V\ — as a natural pro- 
jection of the solution y^ g (t, ft) of the Hamilton-Ehrenfest system (|3.2p onto the phase space, i.e. we set 

Z(t,h) = Z 2 (t,h) = (P 2 (t,h),X 2 (t,h)). Define the function S(t,h) (second "parameter" of the class V%) as 
an analog of the classical action along this (z = Z 2 (t, H),t £ [0, T]) phase trajectory by a standard formula in 
a classical Hamiltonian corresponding not to the main, but to the total symbol [16] 7i 2 t (z, t) of the quantum 
Hamiltonian in (|4.4p . which, in view of the estimates A 2 (t, ft) = 0(ftJ"l/ 2 ), has the form 



Hl(z, t) = H(z, t) + kV(z, w, t) 
As a result, we have 



w=Z{t,<£ T ) 2 



V ww (z,w,t) 



w=Z{t,€ T ) 



A 2 (t,Cr) 



(4.7) 



s(t,n) = s^(t,h) = {(pW(T,h),xW(T,h)) -nW(zW( r ,h),T)}dT. 



(4.8) 



Introduce the notation 



S (2) (x, t, ft) = 5 (2) (t, ft) + (PW ( t , ft), x - (t, ft)) . 



01' 



Expand the operators Ti(z,t) and — — V(z,w,t)\ w=Z (2)/ t h \ in Taylor series of order two in powers of the 

g w a „ 

operator Az = z — Z^ 2 '(t,H) with remainder terms PJ^ and R\ , respectively, and represent the operator 
—ihd/dt in the form 



-ihd t = A + B, 



4 thereby stressing the dependence of the solution of equation 1 14. 31 1 (via the coefficients of the operator ZA 2 )(t, (£)) on the 
solution g(t, h, £) of the Hamilton-Ehrenfest equation 113.21 1 



B = -(P (2) (t,h),X (2 \t,h)) +sW(t,h)- (Z ( V(t,h), JAz), (4.9) 

A = -ihd t - s (2) (t, h) + (t, h), (t, h)) + (z& (t, n), JAz) = d(h). 

Substituting the resulting expressions in equation (|4.3[) . in view of the estimates in % — ► 0, for the operators 
-R3 and -R3 applied to functions of the class V\, we obtain 

{-tfi^ + ho(t, £t)}<s> = 0(h 3 / 2 ), (4.10) 

where 

ho(t, C r )) = -^ (2) (i, ft) + (^ (2) (i, ft), ^ (2) (i, ft)) + (^ (2) (*, ft), ^Az) + I(Az, C r ))A5), (4.11) 

55«(t,Cr)) = [W»(2,*) + >fV r w(« ) «',t)] • (4-12) 

z=to=Z( 2 >(t,ft) 

This is a Schrodinger equation with the Hamiltonian quadratic in the operators p and x. 

5 Trajectory- coherent states 

of the nonlocal Gross— Pitaevsky equation 

The solution of the linear Schrodinger equation with a quadratic Hamiltonian is well known. For our 
purposes, it is convenient to take for the basis of solutions of equation (|4.10[) the semiclassical trajectory- 
coherent states (TCS's) of this equation [23,24]. By Statement [T] these states are asymptotic (mod ft 3 / 2 ) 
solutions of the problem (|2.14p - (|2.12p if the function tp (x,H) (|2.12[) coincides with the TCS at the time 
zero. We shall also call these solutions trajectory- coherent states of the nonlocal Gross-Pitaevsky equation 
(|l.ip - (|1.3p . We now give their explicit form and some properties, which will be used below to solve the 
problem (|2.14p - (|2.12p with an arbitrary initial condition of the class of functions V\. 

Let us use the "momentum" and "coordinate" components of the solution of the Floquet problem (|3 . 1 7|) , 
(I3.25P to compose n x n matrices 

a(t,£ T )) = (W(t,£ T )),Z(t,£ T ))y, 
B(t) = {W x {t), W 2 (t), W n (t)), C(t) = {Z x {t),Z 2 {t), . . . , Z n {t)). (5.1) 

The linear Hamilton system (|3.17p is known [25] to conserve the standard simplectic structure to — dpAdx 
of the phase space R 2 ™ x , and, hence, the skew-scalar product (Ja, b) = u(a, b) of any two complex solutions 
a(t), b(t) 6 C^ z of system p. 171) does not depend on time. (Here, C^ z denotes the complcxification R 2n x 
with complex coordinates w G C™, z G C™, a = (w, z).) From this reasoning, in view of (|3.18p . we obtain 

CT(t)B(t) - B T C(t) = \\{ ai (t), aj (t)}\\ nxn = 0, (5.2) 

j.(C+B(t) - B+C{t)) = ^||{oi(*),o5(t)}|| nXB = I, (5.3) 

Conventional reasoning (see, e.g., [15]), in view of (|5.3p . lead to the statement that the matrix C(t) is 
nondegenerate, det C(t) ^ 0, and the imaginary part of the matrix 

Q(t) =B(t)C- 1 (t) (5.4) 

is positive for t G [0,T], T > 0: lmQ(t) > 0. Moreover, by virtue of (15. 2p . the matrix Q(t) is symmetric, 
QT(t) = Q(t). 

Let us now fix the continuous branch of the root of detC(i), t G [0, T], assuming, for instance, that 
Arg,/det C(0) = 0, and dete rmine the function 

<4 2) (f,i,ft) = \Q,t,<Z T ) = -^B=exp\l:\s^(t,h) + (P< 2 \t,h),Ax) + hAx,Q(t)Ax)}\, (5.5) 

s/ det C(t) ^iii z J J 

where Ax = x - (t, h), {P^ (t, ft), X^ (t, h)) = Z®\t,h)\ t G [0,T] is the phase trajectory by virtue of 
the second order Hamilton-Ehrenfest system (|3.17p ; S^ 2 \t,t£) is defined in (|4.7p . and Nq(H) = (7rfi) - "/ 4 is a 
normalizing constant: (t,0|0,i) = 1. 

Theorem 3 The function $q (0,t,h) (|5.5|) is an (exact) solution of the Cauchy problem, for equation (|4.10p 
with an initial condition of the form 

$[, 2) (f,t,ft)| t= o = No{h)e^ { P ' S ^ a)+ i {s ~ So ' B{0)c ' lms - s ° m , (5.6) 



where (po,Xo) G R 2 ,™,,; the (n x n) complex matrices B(t) and C{t) are defined in (|5.ip . 



Proof. We seek a solution of the linear equation (|4.1ip in the class of functions T>t(S( 2 >(t, h), Z^'{t, H)) in 
the form of a Gaussian packet 



(x, t, h) = exp { l - [S^ (t, h) + (PW {t , h), Ax)] I exp U < Af '^ Af > U (t ), (5.7) 



where the complex (n x n) matrix Q(t) is symmetric and Im Q(t) > 0. Substituting (|5.7[) in (|4.10p and 
equating the coefficients of the powers of the operator Ax k , k = 0, 2, to zero, we obtain a linear equation in 
the function tp(t) and a (matrix) Riccati equation in the matrix Q(t), respectively: 

^ + ^Sp[f) px (t, C T )) + £ pp (i, Cr))Q(*)]v = 0, (5.8) 

Q + «xx(t, Cr) + Q{m P x{t, Cr) + fixp(t, Cr)Q(t) + Q{t)fi PP (t, € T )Q(t) = 0. (5.9) 

The ordinary change of variables Q(t) = B(t)C _1 (i) (see, e.g., [15]), provided that Q(Q) = B^C^O), 
Reduces the problem of constructing the desired complex solution of the equation (|5.9[) to the problem 

* ) - ^ (?) ~ n zt^jtt^ 

which, by virtue of l|5.1[) . is reduced to system (|3.1T[) . By virtue of the second equation of (|5.10[) . we have 

c = [z px (t, ©r)) + fl w (t, e: T ))Q(*)]c, 

where Q(t) is a solution of (|5.9p . Thus, by virtue of Liouville's lemma, we find 

t 

detC(i) = exp J Sp[9) px {t, € T )) + <o pp (t, € t ))Q(t)]cIt, 
o 

and, hence, by equation (|5.8[) . we have ip{t) = (dot C(t)) -1 / 2 . 

Let us now construct the Fock basis of solutions for equation (|4.10p . To do this, we find for this equation 
the symmetry operators a(t, £t)), linear in operators Az, in the form a(t, £t)) = N a (b(t, Cr))> Az), where 
jV a is a constant, 6(t) = b(t, £t)) is the complex 2n vector to be determined. By the equation 

-is^P + [flo(t,CT)),a(*)] = o, 

which determines the operator a(t), in view of the explicit form of 9jo(t, <£t)) (|4.11|) . we obtain 

-ih(b(t),Az) +ih(b(t),zM(t,h)) + 
+ [{-S^(t, h) + (pW(t, h),i^(t, h)) + (ZW(t, h), JAz) + ~(Az,$j zz (t, £t))Az)}, (b(t),Az)] = 0. 



From this equation, by virtue of the commutation relations [Azj, Azk] — ifrJjk, j,k — l,2n, we find 
—ih(b(t),Az) + ih(Az,S) zz (t)Jb(t)) — 0, and, hence, b = Sj zz (t,€T))Jb. Using the notation b(t) = —Ja(t), 
we obtain for the vector a(t) the system in variations (|3.17p . Thus, the operator 

a(t) = a(t, €t)) = N a (b(t),Az) = N a {a(t), JAz) (5.11) 

is the symmetry operator for equation (|4.10[) if the vector a(t) — a(t, Ct)) is a solution of the system in 
variations (|3.17j) . Let a(t) and b(t) be symmetry operators corresponding to two solutions of the system in 
variations, a(t) and b(t), respectively. It is easy to check that 

[6(t), S(t)] = ihN a N b {a{t),b(t)} = ihN a N b {a{0),b(0)}, (5.12) 

the last equation being a corollary of the Hamiltonian character of system (I3.17[) . 

By formula (|5. 1 lj) . associate the vectors a*(t) with the "creation" operators a^(t) and the vectors Oj(i) 
with the "annihilation" operators %•(£), setting N a = {2h)^ 1 / 2 . Then, by virtue of formulas (|5.12p .for the 
operators a~^(t), a,j(t), j = l,n, the following commutation relations, canonical for boson operators, are 
valid: 

[a j (t),a k {t)] = [af(t),a+{t)]=Q, [&j(t),a+{t)] = S jk , j,k = T^i. (5.13) 



(2) 

Statement 2 The function % '(x,t,H) = \0,t,£ T )) iTOl ) is a "vacuum" trajectory-coherent state: 

a J (t)|0,t,Cr))=0, J = M. (5.14) 

Proof. Applying the annihilation operator a,j(t) to the function |0, t), we obtain that fij(0, t)|0, t) = 
|0, t)[(Zj(t), Q(t) Ax) — (Wj(t), Ax)]. From this (|5.14[) immediately follows since, by definition and in view 
of the properties of the matrix Q(t), we have Q(t)Zj(t) = B(t)C~ 1 (t)Z j (t) = Wj(t). 

Let us now define a countable set of states \v, t, £x)) (exact solutions of equation l|4.10p ) as the result of 
the action of the birth operators on the vacuum state |0,t, £t)) (I5.5[) : 



(5.15) 



1 n 1 

$( 2) (z,i,fi) = \v,t,€r)) = -AcL+(t,£ T ))y\0,t,<L T )) = TT — (a+(t, |Q, t, € T )) 



k=l 



( 2") 

The functions '(x,t,H), v € 7/f constitute the Fock basis of solutions of the (linear) equation (j4. 1 0|) : 
by formulas (|5.14p . (|5. 13[) . standard calculations check whether this set of functions is orthonormalized, 
= <W' f v,v' <E Z?., and the proof of its completeness follows, for example, from the results 
presented in [26]. Thus, from this reasoning and Statement [T] we arrive at the following theorem. 

Theorem 4 Let the symbols of the operators Ti.(t) and V(t, in (|1.1[) - (|1.3[) satisfy the conditions of As- 
sumption [7] and the conditions of Theorem be fulfilled. Then, for any v 6 U\ the function 



y u (x,t,K) = <f> u (x,t,h,<£%.) 



(5.16) 



is an asymptotic (to within 0(fi 3//2 ) ; h — > 0) solution of the nonlocal Gross-Pitaevsky equation (| 1 . lj) — (| 1 . 3|) 
with a quasiperiodicity condition (II. 6p . where 



(5.17) 



4 2 >(mod U ) = ~s(r ) fi) + s^n k (i 1t + .), 

^=1 



TVie constants are determined by the equation 



where 



g(t, C T ) = fl o(^o) + 0(^ 3/2 ) = ™<^)|fl|^) + 0(h?' 2 ) 

t=s \\ljj\\ 2 



Mx,n) = ^\x,0,h), 



(5.18) 
(5.19) 



and <£>!, (x,t,fi) are defined in (|5.15j) . 



By the periodicity condition (|3.24[) for the solutions of the Hamilton-Ehrcnfest system and quasiperiod- 
icity condition (|3 . 25[) for the solutions of the system in variations it follows that 



a k (t + T) = e m * T a k {t), a(t) 
>(t + T) 



exp 



2h 

n 



(a(t),JAz) k = l,n, 
a+(t)Y, 



fc=l 



detC*(t + T) = exp [i^TQ fc i/ fe ] detC(t), Q{t + T) = Q(t). 



(5.20) 
(5.21) 

(5.22) 



fc=i 



From the relations 



|0,t) = |0,t,Cr)) = 



{^[S(t, K) + (P(t, h), Ax) + \(Ax, Q(t)Ax)] }, 
t 

S(t,H)= I dt{{P(t,H)X(t,K)) -H(z,t) - xV ww (z,w,t) 



x exp 



Vw 



,(z,w,t)A 2 {t,<tT)]} 



Z=w = Z(t.£rr). 



(5.23) 



(5.24) 



which determine the vacuum trajectory-coherent state, the definition of semiclassical TCS's [5] 

1 " 1 

$ v {x,t,h,£ T ) = M) = W,tXr)) = ^=(a+(t,£TW\0,tX T )) = J] ~r=\ ^ (*> \0,t, C T )). 



fc=l ^ 

(5.25) 

and the quasi-periodicity conditions (|5.21[) it follows that 

M + T) =e~ i£ ^ T ' h \v,t), (5.26) 

where is defined in (|5.17p . Hence, the functions ^ v (x,t,h) (|5.16| also satisfy the quasi-periodicity 
condition (|5.26p and constitute the quasi-energy spectral series [^ u (x, t, h), £„ ] of equation (jTTTJ) which 
correspond to the phase curve z = Z (t). 

6 Semiclassically concentrated solutions of the nonlocal Gross— 
Pitaevsky equation (the principal term of the asymptotics) 

The asymptotic solutions of the Cauchy problem (|2.14[) - (|2.12p (where ipo has the form of (|5.19p ) that 
have been constructed in the previous section are a special case of the semiclassically concentrated ( mod H 3 ^ 2 ) 
solutions of equation (|l.ip - (ll.3p . In case of arbitrary initial conditions ipo(x, H) (12. 12|) belonging to the class 
Vft{Z^(t,H),S^\t,K)\, t = s, the principal term of the asymptotics \I/q 2 (x, t, h) of the nonlocal Gross- 
Pitaevsky equation (|1.1|) - (|1.3|) is determined by expanding the solutions of equation (|4.10p in a series over 
the Fock basis \u,t, <L T )) ([5TT5|) . v G Z": 

OO 

^ 2) (f,t,fc) = aiM.Mifo)). c v = (_ T (^o),M^o>. 

|v|=0 

Constants Cr(V'o) are determined by the equation (|3.1ip . For the subsequent calculations, it is convenient 

(2) 

to represent the solution -0q (x, t, K) in terms of the convolution of the initial condition ipo with the Green's 

( 2^ — » 

function G {£, y*, t, s, Crj) of equation (|4.10p . The Green's function for quadratic quantum systems is 
well known (see, e.g., [13,27]). For completeness we give its explicit form in a convenient representation 
allowing us to clearly demonstrate the nontrivial dependence of the evolution operator of the associated 
linear equation (|4.10p on the initial conditions for the input nonlocal Gross-Pitaevsky equation. 

/ 2~j _____ 

By the definition of Green's function G , for equation (|4.10p we have 

[-ihdt + ft (t, £t))]G^ 2) (a* y, t, s, € T )) = 0, < s < t; 

(2) (6- 1 ) 
limG '(S,y,t,s,&T)) = 6(£-ff), 

t—>s 

where the operator iOo is defined in (|4.1ip . Denote by Xk(t, (_r), k = 1,2,3,4, the n x n matrices that are 
blocks of the fundamental matrix of the system in variations (15. 10[) : 



<5( 



Let the following conditions be fulfilled: 

detfi pp (t,£ T )) j_0, detX 3 (t-s,€ T ) 7^0, s,t € [0,T], (6.3) 

( 2) 

Then the Green's function Gq '(x^y^t, s, £t)) has the form 

+(pW(t, h),Ax) - (po, (y- 5b)) - \{{y~ x ),Xi(At, £ T )A^(At, € T )(y- x ))+ 

+{AS, X^(At, <L T ){y- x )) - i(Aa;, X^(At, £ T )X 4 {At, £ T )A£)] }. (6.4) 



Here, At = t — s, Z^(t,H) = (p( 2 \t, H), X^ 2 \t, H)) is the projection of the solution of the second order 
Hamilton-Ehrcnfest system onto M. 2n x . 



Remark 1 If condition (|6.3[) is not fulfilled, the solution of the problem (|6.1[) has a somewhat different form 
than (|nill> (see, e.g., [27]). 

From this reasoning and Statement [1] we arrive at the following theorem. 

Theorem 5 Let the symbols of the operators Ti.(t) and V(t,^) in (|1.1[) - (|1.3[) satisfy the conditions of As- 
sumption [7] and let the conditions of Theorem [3| and conditions (|6 . 3[) be fulfilled. Then the function 

4 2) (x,t,h)^U^(t,0)^ , t€[0,T], (6.5) 

where Uq (t,0) is the evolution operator of the zero order associated Schrddinger equation 1)4. 10p with a 
kernel G 2 (x, y, t, 0, €(ipo))) (|6 .4(1 . is an asymptotic (to within 0(h 3 ^ 2 ), h — > 0) solution of the nonlocal 
Gross-Pitaevsky equation - (|1.3p with the initial condition 

*|«=o = M 2 )- 

Corollary 5.1 Operator Uq (T) = Uq (T, 0) zs i/ie monodromy operator the Gross-Pitaevsky equation 

_______ "^"(2) 

(|1.1|) (|1.3|) . Functions ^ u (x,t,K) (I5.16P are eigenf unctions of operator Uq (T): 

U {2) (T)*„(f , t, K) = e-^^/^^f , t, fi). 



7 Geometric phases of trajectory-coherent states 

We now turn to the calculation of the Aharonov-Anandan phase corresponding to the quasi-energy states 
(I5.17p . To do this, we use formula (| 1 . 1 3[) by which, neglecting the terms of order 0(H), we obtain the following 
expression of the desired phase if the phase curve is periodic: 



Is. 



£ { U 2) T 1 



T 

[ dt\n{t)+kV(t) + hp(sj zz {t)A 2 ) + JSp[VL, w (Z(i,R),w,i)A 2 ]}. (7.1 

J ^ Z V / Z u}—Z{t,h) J 



h h 

Using (|3.22p , we can readily show that 
Sp([H zz (t) + xV zz (t)]A 2 (t) 



Sp Re 



C(t)D^B + (t)-B(t)D^C + (t) 



n 

HReJ2U + -){a k (t),am}, 



k=l 



(7.2) 



where, in our case, D v 1 = diag(2^i + 1, . . . , 2v n + 1). Note that the quantity under the summation sign on 
the right of (|7.2p . by virtue of (|3 . 1 8|) . is real. Hence, the sign Re can be omitted. Introduce the notation 
ao(i) = Z (t), where Z(t, h) = Z (t) + %Z\{t) + 0(/i 3/2 ) are defined in ([3TT5]) . Then, substituting the explicit 
expression for quasi-energies £„ (|5.17p in (|7.ip and using relation (|7.2p . we find 

T T T 



= \ I dt(P(t),X(t))- I A{Z (t),Zi(t)}-X; 



7£„ 



fc=l 



1 

Vk+ 2 



TO fe + i / ,//{«,,(/ I. r/ A (/l} 



If now we introduce, instead of the Floquet solutions, the T-periodic vector functions 

a k (t) = e-' inkt a k {t), ~a k {t + T) = a k (t), 
then formula (17.311 takes the form 



(7.3) 



(7.4) 



l£„ 



1 



dt{P(t),X(t))- / d*{ao(t),Zi(t)} ^ 



1 



fc=i 



dt{a k (t),a* k (t)}. 



Using (|3.23p , we obtain the following expression for the Aharonov-Anandan phase (|7.5[) : 



(7.5) 



l£ v 



=- 1 dt(P(t),x(t))--Y: 



fc=i 



^ fe+ 2 



dt{a k ,a* k }- 



-Re |^^{a (t),oKt)}£(^ + i) | dr{o?(r), d z }{a* k (r), j[H zz (z, r)- 
o ' =1 * =1 o 

+ xV zz (z,w,t) + kV ww (z,w,T)]a k {T)\ +b k {T)e mk 



(7.6) 



Thus, the the Aharonov-Anandan phase j£ v corresponding to the quasi-energy TCS's (|5 . 16[) is com- 
pletely determined by two geometric objects: the closed phase trajectory Zo(t) of the Hamilton- Ehrcnfcst 
system, stable in the linear approximation, and the complex germ r n (Zo(t)) composed of n linearly inde- 
pendent Floquet solutions of the system in variations (|3.17[) . 

Now consider Eqs. (jl.lj) — (11.3|) in a 3-dimensional space with operators Tt(t), V(t, *S>) of the form 

H(z,t) = ^-f - e(E(t),x) + 
2m 2 

{x-y) 2 - 



V(z, w, t) = V(x — y) = Vq exp 



2 7 2 



(7.7) 
(7.8) 



The external field in the linear operator (|7.7[) is electric field E(t) — (E cos cut, E sin uit, 0) periodic in time 
with frequency u>, and the field of an isotropic oscillator with potential kx 2 /2, k > 0. 
The periodic solution (|3.16|) corresponding to (|7.7p . (|7.8p accurate to 0(H 3 / 2 ) 



Z (t,H) = (P (t,h),x (t,h)), 

Pa(t, h) = ( — mw( smut, mu( cos tot, 0^ , 
Ao (t , h) = (£ cos tot , £ sin cot , 0) . 

Here £ = eE/[m(oj 2 ~ to 2 )]. 

The Floquet solutions of system (]3 . 1 7[) corresponding to (|7.7p . (|T.8[) accurate to 0(h 3 ^ 2 ): 



(7.9) 



Ol(t) 

Oa(*) 
a 3 (t) 



st ( 5s ,o,o,--,o,o) T , 

■*(0,ff B) 0,0,--,0)T ) 
v 5s 

st (o,o, 5s , 0,0,--) T . 

v 5s 



(7.10) 



Here 



5s = V m ^s; Wnl = 



The spectrum of quasi-energies £u (|5.17p is given by the relation 



IxVbl 

7717 2 



OJ 



h = k\\^\\ ; 77 = sign (i<Vb). 



?(2) 



2w s 



)( 



v\ + - 



2uj s 



1 

f 2 + ~ 

2 



2uj s 



U 3 + 77 



0(»' 



3/2x 



The Aharonov-Anandan geometric phase (|7.6p is 

Tu 2 



me + 0(h 1/2 ). 



(7.11) 



(7.12) 



We have found the explicit expressions for the Aharonov-Anandan phase in semiclassical approximation 
accurate to 0{H 1 / 2 ) for the quasi-energy states which are governed by the GPE and belong to the class of 
trajectory concentrated functions. It is of interest to consider the expressions (|7.6p in the adiabatic limit 
T — > 00 and to observe in what sense the Aharonov-Anandan phase transforms to the Berry phase in 
the background of formalism of the GPE. The fact is that the splitting of the total phase, gained by a 
wave function, into dynamic and geometric parts is carried out by analogy with the case of linear quantum 
mechanics. Then it is not obvious that the expressions obtained for the nonlinear GPE are determined only 
by geometry of the system and do not contain a dynamic contribution due to the nonlinear term in the 
equation. 

In the limit T — > 00 we can neglect the dependence of operators Ti.(t) and V(t, 1 J r ) of the form (jl.ip on 
time t. Therefore, the quasi-energy spectral series (|5.16l) grade into discrete spectral series of the nonlinear 
problem 

[H + xV{^)]9 = 

in the limit T — > 00. These spectral series [28] are localized in a neighborhood of stationary solutions of the 
Hamilton-Ehrcnfcst system (|3.2p , (|3. 16|) . The geometrical phases (|7.12p do not contain quantum corrections 
accurate to O^h 1 / 2 ) in the considered example and do not depend on nonlinear potential. Evidently, it is 
related with the fact that the nonlinearity and the quantum corrections with the given accuracy in H do not 
change solutions (|7.9p of the Hamilton system (|3.16p and, as a result, the geometry of the system. 
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